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Pop Quiz 

i 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0 


y=x 

y=0.6sin(1.7x+l) 



0 
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Pop Quiz 

• Set equations equal 

• x = 0.6sin(1.7x + 1) 
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Pop Quiz 

• Set equations equal 

• x = 0.6sin(1.7* + 1) 

• Solve for x 
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Pop Quiz 

• Set equations equal 

• x = 0.6sin(1.7* + 1) 

• Solve for x 

• No algebraic solution! 
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Numerical Methods 

• Newton's method to the rescue 
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Numerical Methods 


Newton's method to the rescue 
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Fixed Point Iteration (FPI) 

• Given equation of form: 

• x = F(x ) 


• And initial guess x 0 

• Compute iterates x 1 ,x 2 ,x 3 ... 

• *<+i = F(Xi) 
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Fixed Point Iteration (FPI) 

• Equation 

• x = 0.6sin(1.7a: + 1) 

• Let F(x ) = 0.6sin(1.7x + 1) 

• Initial guess * 0 = 0.05 

• Iterate: 

. = F(xo) = F( 0.05) = 0.531 

. x 2 = F(xt) = F( 0.531) = 0.567 




x 5 = f(x 4 ) = F( 0.559) = 0.558 
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Fixed Point Iteration (FPI) 



y=x 

y=0.6sin(1.7x+l) 

FPI 


GDC 


<> 


GAME DEVELOPERS CONFERENCE March 14-18, 2016 - Expo: March 16-18, 2016 #GDC16 


Talk Outline 

• Real world applications of FPI 

• Convergence requirements 

• Power-ups 
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Case 1: Gerstner waves 

• Deep water ocean wave shape [6] 

• Used in Uncharted 3, Disney Infinity: Pirates of the 
Caribbean [1,2] 



GDC 


<> 


GAME DEVELOPERS CONFERENCE March 14-18, 2016 - Expo: March 16-18, 2016 #GDC16 




Case 1: Gerstner waves 


• Problem: Given x, compute y 

• Parametric curve 


• x = t — B sin t 
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Case 1: Gerstner waves 

• Problem: Given x, compute y 

• Parametric curve 

• x = t — B sin t 

• y — A cos t 

• y is a function of t 
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Case 1: Gerstner waves 

• Problem: Given x, compute y 

• Parametric curve 

• x = t — B sin t 

• y — A cos t 

• y is a function of t 

• Need t from first equation 

• Same equation form as pop quiz! 
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Case 1: Gerstner waves 

• Rearrange to get F(t) = t form 

• x = t — B sin t 

• t = x + B sin t 

• Let F(t) =x + #sint 

• Initial guess: t 0 = x 

• Iterate: t t = F(t 0 ), ... 

• Compute y using t n 
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Case 1: Gerstner waves 



Diagram : shadertoy.com/view/MstSRN 
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Case 1: Gerstner waves 



Diagram : shadertoy.com/view/MstSRN 
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Case 1: Gerstner waves 



■D.4 

1.0 O.B 


Diagram : shadertoy.com/view/MstSRN 
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Case 1: Gerstner waves 



t 2 — F(t i) 



-D.H 

M.B 0.8 


Diagram : shadertoy.com/view/MstSRN 
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Case 1: Gerstner waves 



£3 — ^(^2) 



■D.4 

1.0 O.B 


Diagram : shadertoy.com/view/MstSRN 
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Case 1: Gerstner waves 



■D.4 

1.0 OJ 


Diagram : shadertoy.com/view/MstSRN 
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Case 1 Summary 

• FPI can solve hard problems 

• Extremely general 

. t = F{t) 

• More octaves? Sure! 

• Stable - no gradient 
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Case 2: Vector Displacement Maps 


• Simulate the ocean offline instead [1,4] 
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Case 2: Vector Displacement Maps 

• Simulate the ocean offline instead 



[Source: Blender] 
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Case 2: Vector Displacement Maps 



ink 
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Case 2: Vector Displacement Maps 
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Case 2: Vector Displacement Maps 

• Goal: compute height at x 



x 
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Case 2: Vector Displacement Maps 

• Goal: compute height at x 

• Height given at u 



x 


u 
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Case 2: Vector Displacement Maps 

• Need to find u which will be displaced to x 

• Equation: u + D(u) = x 
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Case 2: Vector Displacement Maps 


• FPI to the rescue! 

• Equation : m + D(u) = x 

• Rearrange: u = x - D(u) 

• Let F(u) = x — D(u ) 

• Guess uq = x 

• Iterate: u 1 = x — D(u 0 ), ... 
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Case 2: Vector Displacement Maps 

• Simulate the ocean offline instead 


- 


[Source: Blender] 
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Case 2 Lessons 

• FPI great for searching data 

• Simple to implement 

• Under-iterated result stable 

• Generalises to higher dimensions 


GDC 


<> 


GAME DEVELOPERS CONFERENCE March 14-18, 2016 - Expo: March 16-18, 2016 #GDC16 


Case 3: Distance field terrain 

• Render terrain represented by Distance 
Estimator (DE) 
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Case 3: Distance field terrain 


• Pixel ray 

• p = o + td 

• Want intersection of ray with surface 

• DE(x) = 0 

• Plug in ray: 

• DE(o -h tcT) = 0 

• Solve for t! 

• But this is not our usual form t = F(t ) 
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Case 3: Distance field terrain 

• Easy fix: add t to both sides 

• 0 = DE(o td) 

• t = t - 1 - DE(o -h td) 

• Let = t + DE(o + td) 

• Start at camera: guess t 0 = 0 

• Iterate: t i+1 = t t + DE(o + t t d ) 

• Standard raymarch algorithm! 


GDC 
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A 
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Case 3 Lessons 

• Raymarching is equivalent to FPI 
. Why? 

• FPI is fundamental 

• Any iterative update can be posed as FPI 
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Convergence Requirements 

• Standard form: x = F(x) 

• Converges when |F'(x)| < 1 around solution 

• Study raymarch case: 

• t = t + DE(t ) 

• Converges when |1 + DE'(t)\ < 1 

• Therefore -2 < DE'(t) < 0 




-2 < DE'(t) < 0 






V. 
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Power-ups 

• If flip flops around solution 

• Secant step 

• Jumps over/misses solution? Scale down step. 

• Very common with distance estimators 

• Slow to reach solution? Scale up steps... 

• Over-relaxation for raymarching [5] 
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The Pattern 

• Write down equation with unknown 

• Analytical solve? You're done. 

• Get in form x = F(x ) 

• Iterate x i+1 = F(xi) to taste 
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Summary 

• FPI arises naturally 

• Simple and easy to try 

• Coherent, stable, good for searching data 

• Unlike Newton! 

• If diverges 

• Try to reduce step size 

• Try Secant Step 


GAME DEVELOPERS CONFERENCE March 14-18, 2016 - Expo: March 16-18, 2016 #GDC16 



Summary - Higher Dimensions 

• FPI not limited to ID 

• See [3] for application of FPI to image warping 

• Simple to implement, fast 

• Convergence condition for >1D given 

• Behaviour illustrated in the presence of many solutions 

• Technique described to catch all solutions 


lllk 
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Thank you for listening! 

• And thanks to.. 

• Halldor Fannar, Daniel Zimmermann, Studio Gobo 


• MiSC FlOteS! huwbowles.com/fpi-qdc- 2016 / 

• Contact 

• Email: huw.bowles@qmail.com 

• Linkedln: linkedin.com/in/huwbowles 


